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ABSTRACT 

Exoplanet orbital eccentricities offer valuable clues about the history of planetary systems. Eccen- 
tric, Jupiter-sized planets are particularly interesting: they may link the "cold" Jupiters beyond the 
ice line to close-in hot Jupiters, which are unlikely to have formed in situ. To date, eccentricities of 
individual transiting planets primarily come from radial velocity measurements. Kepler has discov- 
ered hundreds of transiting Jupiters spanning a range of periods, but the faintness of the host stars 
precludes radial velocity follow-up of most. Here we demonstrate a Bayesian method of measuring 
an individual planet's eccentricity solely from its transit light curve using prior knowledge of its host 
star's density. We show that eccentric Jupiters are readily identified by their short ingress/egress/total 
transit durations - part of the "photoeccentric" light curve signature of a planet's eccentricity — even 
with long-cadence Kepler photometry and loosely-constrained stellar parameters. A Markov Chain 
Monte Carlo exploration of parameter posteriors naturally marginalizes over the periapse angle and 
automatically accounts for the transit probability. To demonstrate, we use three published transit 



light curves of HD 17156 b to measure an eccentricity of e = 0.71 



-1-0.16 
-0.09' 



in good agreement with the 



discovery value e = 0.67 ± 0.08 based on 33 radial-velocity measurements. We present two additional 
tests using actual Kepler data. In each case the technique proves to be a viable method of measur- 
ing exoplanet eccentricities and their confidence intervals. Finally, we argue that this method is the 
most efficient, effective means of identifying the extremely eccentric, proto hot Jupiters predicted by 
Socrates et al. (2012). 

Subject headings: planetary systems, techniques: photometric 



1. INTRODUCTION 

Many exoplanets have highly eccentric orbits, a trend 
that has been interpreted as a signature of the dynamical 
processes t hat shape the architectures of planetary sys- 
tems (e.g. 'Juric & Trenmin^ [200l iFord fc Rasiol [20081: 
[Nagasawa & Ida 2011). Giant planets on eccentric or- 
bits are of particular interest because they may be relics 
of the same processes that created the enigmatic class 
of planets known as hot Jupiters: planets on very short 
period ( P < 10 days) orbits th at, unlike smaller plan- 
ets (e.g. [Hansen fc Murravll20i^ . could not have formed 
in situ. Hot Jupiters may have smoothly migrated in- 



ward through the 



Goldreich & Tremain 
20051: llda fc LinI I200i 



disk from which they forme d (e.g. 
1 1980; Ward 1997; Alibert^eLal] 
iBromlcv fc Kenvon. , 



2011, ). Al- 
ternatively, the typical hot Jupiter may have been per- 
turbed by anoth er body onto an eccentric orbit (see 
iNaoz et al1l2012r ). with a star-skirting periapse that be- 
came the parking spot for the planet as its orbit circu- 
larized through tidal dissipation, in itiated by one of sev- 
eral perturbation mechanisms (e.g. iWu fc Murravl 120031 : 
IFord fc Rasiol[2006l : ImTfc Lithwickil20'Tlli ! 
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iSocrates et al.l ()2012l ) (hereafter S12) refer to this pro- 
cess as "high eccentricity migration" (HEM). If HEM 
were responsible for hot Jupiters, at any given time 
we would observe hot Jupiters that have undergone full 
tidal circularization, failed hot Jupiters that have tidal 
timescales too long to circularize over the star's lifetime, 
and proto hot Jupiters that are caught in the process of 
tidal circularization. S12 predicted that the Kepler Mis- 
sion should detect several "super-eccentric" proto hot 
Jupiters with eccentr icities in ex c ess of 0.9. This pre- 
diction was tested bv lDong et al.l (|2012| ) on a sample of 
eclipsing binaries in the Kepler field: in an incomplete 
search, they found 14 long-period, highly eccentric bina- 
ries and expect to eventually find a total of 100. 

As a test of planetary architecture theories, we are de- 
voting a series papers to measuring the individual ec- 
centricities of the Kepler Jupiters to either identify or 
rule out the super-eccentric proto hot Jupiters predicted 
by S12. In this first paper, we describe and demon- 
strate our technique for measuring individual eccentrici- 
ties from transit light curves. Measuring the eccentricity 
of a Jupiter-sized planet is also key to un derstandin g its 
tidal history (e.g. iJackson et al.lf2 008: Ha nsenll2010f) an d 
tidal heating fe.g iMardlind 120071; .Jackson et al.ll2008D . 
climate variations fe.g iKataria et al.ll2011[ ). and the ef- 
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feet of the variation in insolation on the h abitability 
(e.g iSpiegel et ahl 120101 [Dressing et "all I2010D of possi- 
ble orbiting rocky e xomoons detectable by Kepler (e.g. 
iKippinget an[2009( ). 

To date, the measurements of eccentricities of individ- 
ual transiting planets have been made through radial ve- 
locity follow-up, except w hen the planet exhibi ts tran- 
sit timing variations (e.g. iNesvornv et al.l [20121 ). How- 
ever, a transit light curve is significantly affected by a 
planet's eccentricity, particularly if the photometry is 
of high quality: we refer to the signature of a planet's 
eccentricity as the "photoeccentric" effect. One aspect 
is the asymmetry between i ngres s and egress shapes 
(jBurke et all l2007t iKippingj I2008D . The eccentricity 
also affects the t iming, duration, aiid exis t ence of sec- 
onda ry eclipses (jKane fc von Brauiil 120091 : iDong et aD 
l2012f ). The most detectable aspect of the photoeccentric 
effect in Kepler photometry for long-period, planet-sized 
companions is the transit event's duration at a given or- 
bital period P, which is the focus of this work. 

Depending on the orientation of the planet's argu- 
ment of periapse (a;), the planet moves faster or slower 
during its transit than if it wer e_on_a circu l ar or- 
bit w ith the same orbital period (Barnes 20071 iBurkd 
20081 Ford. Quinn. a nd Veras 2008, hereafter FQV08; 



Moorhead et al.l 120111) . If the transit ingress and egress 



durations can be constrained, the duration aspect of the 
photoeccentric effect can be distinguished from the effect 
of the planet's impact parameter (6), because although 
6 > shortens the full transit duration (123, during 
which the full disk of the planet is inside the disk of 
the star, i.e. from second to third contact), it lengthens 
the ingress/egress duration. Therefore, with prior knowl- 
edge or assumptions of the stellar parameters, combined 
with measurements from the light curve of the planet's 
period and size (Rp/Ri,), one can identify highly eccen- 
tric planets as those moving at speeds inconsistent with 
a circular orbit as they pass in front of their stars (see 
also §3 of IBarnesl [2007[ §3.1 of FQV08). 

B arnes! ()2007t ) presented the first comprehensive de- 
scription of the effects of orbital eccentricity on a transit 
light curve, including that a short transit duration cor- 
responds to a minimum eccentricity, contingent on the 
measurement of b and of the host star's density. IBurkd 
([^08) discussed the effect of orbital eccentricity on tran- 
sit detection and on the inferred distribution of plan- 
etary eccentricities. FQV08 laid out the framework for 
using photometry to measure both the distribution of ex- 
oplanet eccentricities and, for high signal-to-noise tran- 
sits of stars with known parameters, the eccentricities of 
individual planets. They derived expressions linking the 
orbital eccentricity to the transit duration and presented 
predicted posterior distributions of eccentricity and uj for 
a given ratio of: 1) the measured total transit duration 
(i.e. from first to fourth contact, including ingress and 
egress) T14 to 2) the T14 expected for a planet on a circu- 
lar orbit with the same &, stellar density p^, and P. Then 
they showed how the distribution of planetary transit du- 
rations reveals the underlying eccentricity distribution. 
FQV08 focused on the possibility of measuring the ec- 
centricity distribution of terrestrial planets, which has 
implications for habitability. Here we will show that the 
technique they describe for measuring individual planet 
eccentricities is particularly well-suited for Jupiter-sized 



planets. 

The work of FQV08 was the basis for several recent 
analyses of high-precision light curves from the Kepler 
mission that have revealed information about the eccen- 
tricity distribution of extra-solar planets and the eccen- 
tricities of planets in multi-transiting systems. By com- 
paring the distribution of observed transit durations to 
the distributi on derived from model p opulations of eccen- 
tric planets, iMoorhead et al.l (|201ll ) ruled out extreme 
eccentricity distributions. They also identified individual 
planets with transit durations too long to be consistent 
with a circular orbit; these planets are either on eccentric 
orbits (transiting near apoapse) or orbit host stars whose 
st ellar radii a. r e sign ificantly underestimated. 

iKane et al.l (|2012D used the distribution of transit du- 
rations to determine that the eccentricity distribution of 
Kepler planets matches that of planets detected by the 
RV method and to discover a trend tha t small planets 
have less eccentric orbits. In contrast, iPlavchan et al.l 
()2012|) found that the distribution of eccentricities in- 
ferred from the transit durations is not in agreement with 
the eccentricity distribution of the RV sample; they sug- 
gested that the difference m ay be due to e r rors i n the 
stellar parameters. Finally, iKipping et al.l ()2012[ ) pre- 
sented a method that they refer to as Multibody Astero- 
density Profiling to constrain eccentricities of planets in 
systems in which multiple planets transit. They noted 
that one can also apply the technique to single transit- 
ing planets, but discouraged doing so, except for planets 
whose host star densities have been tightly constrained 
(e.g. by asteroseismology). FQV08 recommend measur- 
ing eccentricities photometrically only for planets with 
"well-measured stellar properties" but also point out the 
weak dependence of eccentricity on stellar density. 

In this work we apply the idea first proposed by FQV08 
to real data and demonstrate that we can measure the 
eccentricity of an individual transiting planet from its 
transit light curve. We show that this technique is par- 
ticularly well-suited for our goal of identifying highly ec- 
centric, giant planets. In §2, we show that even a loose 
prior on the stellar density allows for a strong constraint 
on the planet's orbital eccentricity. In §3, we argue that 
Markov Chain Monte Carlo (MCMC) exploration of the 
parameter posteriors naturally marginalizes over the pe- 
riapse angle and automatically accounts for the transit 
probability. We include both a mathematical and practi- 
cal framework for transforming the data and prior infor- 
mation into an eccentricity posterior. In §4, we measure 
the eccentricity of HD 17156 b from ground-based tran- 
sit light curves alone, finding good agreement with the 
nominal value from RV measurements. We also measure 
the eccentricity of a transit signal injected into both short 
and long cadence Kepler data and of Kepler Object of In- 
terest (KOI) 686.01 from long-cadence, publicly-available 
Kepler data, finding an eccentricity of e = 0.621q :[4. In 
§5, we present our program of "distilling" highly eccen- 
tric Jupiters from the KOI sample and we conclude (§6) 
with prospects for further applications of the photoec- 
centric effect. 

2. PRECISE ECCENTRICITIES FROM LOOSE 
CONSTRAINTS ON STELLAR DENSITY 

To first order, a transiting planet's eccentricity and its 
host star's density depend degenerately on transit light 
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curve observables. IKipping et al.l (I2012D harnessed the pow er of muhiple planets transiting the same host star to 
break this degeneracy fsee also iRagozzine fc Holniai]|2010f) . Yet, as FQV08 first pointed out, although the transit 
observables depend on the stellar density, this dependence is weak (the ratio of the planet's semi-major axis to the 

stellar radius a/Ri, ex pj ). Thus a loose prior on the stellar density should allow for a strong constraint on the 
eccentricity. 

In the limit of a constant star-pla net distance duri ng transit and a non-grazing trans it (such that th e transit is 
approximately centered at conjunction). iKippingl (|2010al ) derived the following expression (|Kipping|l2010al Equations 
30 and 31) for T14, the duration from first to fourth contact (i.e. the total transit duration including ingress and 
egress), and for T23, the duration from first to third contact (i.e. the full transit duration during which the full disk 
of the planet is inside the disk of the star) : 



Ti 
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(1) 



where P is the orbital period; e is the eccentricity; lo is the argument of periapse; R^, is the stellar radius; S — (Rp/Ri,)'^ 
is the fractional transit depth with Rp the planetary radius; a is the semi-major axis; and i is the inclination. By 
combining T14 and T23, we can rewrite Equation ([T]) as 
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Using the small angle approximation, which is also used bv IKippm^ (|2010aD . allows us to group the transit light 
curve observables on the right-hand side: 
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where 
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The g notation is inspired bv IKipping (|2010aD and IKipping et all ()2012fl 's variable ^, for which "if = g'^. Dynamically, 
g is the ratio of the planet's velocity during transit (approximated as being constant throughout the transit) to the 
speed expected of a planet with the same period but e = 0. Note that w is the angle of the periapse from the sky 
plane, such that w = 90° corresponds to a transit at periapse and uj = —90° to a transit at apoapse. For a given 
P and 5, T14 and T23 are shortest (longest) and g largest (smallest) when the planet transits at periapse (apoapse). 
Moreover, if we approximate sini = 1, we can rewrite Equation ([3]) as: 
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Finally, using Kepler's third law and assuming that the planet mass is much less than the stellar mass {Mp ^ Mi,) 
the transit observables can be expressed in terms of the stellar density p^,: 



P*(e,a;) = g{e,uj) ^pc 



where 



Pcirc = p*(e = 0) 



251/4 



V ^1^4 ~ ^2^3 



3P 



(6) 
(7) 



Although Equation |6] was derived under several stated approximations, the relationships among p*, e, and w are key 
to understanding how and to what extent we can constrain a transiting planet's eccentricity using a full light curve 
model. Because g{e,uj) is raised to such a large power, a small range of g[e,u>) corresponds to a large range in the 
ratio p*/pcii-c, i-e. the ratio of the true stellar density to the density measured from fitting a circular transit light 
curve model. For instance, the assumed value of p* would need to be in error by two orders of magnitude to produce 
the same effect as a planet with e = 0.9 and to — 90°. Thus the pcirc derived from the transit light curve strongly 

constrains g, even with a weak prior on p*, because g (x pi. 
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2.1. Constraints on pcirc from the light curve: common 

concerns 

One might worry that long-cadence data, such as the 
30-minute binning of most Kepler hght curves, cannot 
resolve the ingress and egress times sufficiently to con- 
strain a/Ri,, or equivalently Pcirc- In other words, one 
might worry that a/ Ri, is completely degenerate with b, 
and hence that the denominator of Equation ^ is un- 
constrained. This is often the case for small planets. 
However, Jupiter-sized planets have high signal-to-noise 
transits and longer ingress and egress durations (due to 
the large size of the planet). See §2.1 of FQV08 for an 
analysis of how the precision of Kepler data affects con- 
straints on the total, ingress, and egress durations. 

Furthermore, even if the ingress is unresolved or poorly 
resolved, it is often impossible for the impact parameter 
b to account for the short duration of a highly eccen- 
tric, Jupiter-sized planet's non-grazing transit. The max- 
imum non-grazing impact parameter is 1 — Rp / R^, < 0.9 
for a Jupiter around a Sun-like star. Imagine that 
an eccentric planet transits at zero impact parame- 
ter (i.e. travels across 2Rp -\- 2i?*) at speed g. If 
we instead assume that planet is transiting at its cir- 
cular speed g = 1 across the short chord of length 

(2y(^7T"RjO^~~(^iIJ^ir^J^^^^^^^iAF)i the required im- 
pact parameter would be: 



^largc enough 



(i + 5i/2)yT^T7P 



(8) 



For g = 2.38 (corresponding to e = 0.7, a; — 90°) and 
,51/2 ^ Rp/R^ = 0.1, b would need to be « 0.998, 
which would be inconsistent with a non-grazing transit. 
In contrast, a planet with Rp/R^ = 0.01 would have 
^largc enough ~ 0.917, Consistent with the b < 0.99 neces- 
sary for a non-grazing transit. We note this effect simply 
to hig hlight a constraint that arises naturally when fit- 
ting a iMandel fc Agoll (|2002f) transit model to a hght 
curve. 

Additionally, with a properly binned model (as dis- 
cussed in Kipping 2010b, who advocates resampling the 
data times, computing a model light curve, and then 
smoothing to match the data cadence), multiple tran- 
sits allow for constraints on the ingress and egress, 
even if they are poorly resolved in a single transit. 
We demonstrate eccentricity measurements using long- 
cadence data in §4.2. 

Another concern regards the degeneracy of a/Ri, and 
b with the limb-darkening parameters. Limb darken- 
ing causes the shape of the transit to be rounded in- 
stead of flat, potentially causing confusion between the 
full transit and the ingress/egress. However, in practice 
we find that it makes little difference whether we freely 
vary the limb darkening parameters or impose a normal 
prior based on the stellar parameters (e.g. the coeffi- 
cients computed for the Kepler bandpass bv lSindl2010D . 
FQV08 also find that limb darkening does not have a sig- 
nificant effect on the other parameters, as demonstrated 
through tests on simulated light curves (see FQV08 §2.1 
and FQV08 Figure 5). 

Finally, one might worry about dilution by light from a 
near by or background st ar blended with the target star 
(see iJohnson et ahl 120111 for a Kepler example). Dilu- 
tion would cause Rp/R^, to appear too small. Consider 
the impact that dilution would have on the derived pa- 



rameters of an eccentric planet transiting near periapse. 
The ingress and egress durations would be longer than 
expected, and the inferred maximum impact parameter 
to avoid a grazing orbit (i.e. 1 — Rp/ R^) would be too 
large. Both of these effects would caused the planet's or- 
bit to appear less e ccentric (or, equivalently , for pcirc to 
appear smaller; see lKipping fc Tinettill2010l for a formal 
derivation of the effect of blending on the measurement 
of a/Ri,). Therefore, dilution would not cause us to over- 
estimate a planet's eccentricity, if the transit duration is 
shorter that circular. Moreover, because pcirc depends 
only weakly on the transit depth (Equation [7|), the effect 
of blending on the eccentricity measurement is small. We 
quantify this effect through an example in the next sub- 
section. 

Furthermore, if we were to mistakenly attribute an 
apparently overly-long transit caused by blending to a 
planet transiting near apoapse, the resulting false eccen- 
tricity would be quite small. Imagine that the planet is 
on a circular orbit, but that the blend causes us to mea- 
sure pcirc = (1 ^ f)p*i where < / <C 1. The inferred g 
would he g = [pcac/ P*Y^^ « 1 — //3, very close to the 
true g = 1 of the circular orbit. 

2.2. Constraints on eccentricity 

From Equation [51 it might appear that e and lo are in- 
extricably degenerate for a single transiting planet. Cer- 
tainly, if pcirc is consistent with p^, any eccentricity is 
consistent with the transit observables. However, a nom- 
inal value of Pi, smaller than pcirc translates to a mini- 
mum eccentricity eminj the value obtained by as suming 
the p lanet transits at periapse (w = 90°; see also iBarned 
120071 §3; Kane ct al. 2012, §4). Conversely, a value of 
Pi larger than pcirc corresponds to an Cmin obtained by 
assuming the planet transits at apoapse (w — —90°). 
Therefore, we can easily identify planets with large eccen- 
tricities. A full MCMC exploration provides a confidence 
interval that shrinks as e — > 1, as we discuss in detail in 
§3. For example, consider a planet with an eccentricity of 
0.9 that transits at semilatus rectum (w = 0). Based on 
the transit light curve observables, we would deduce that 
it has an eccentricity of at least emin = 0.68. A planet 
transiting at semilatus rectum with e = 0.98 would have 
a deduced Cmm = 0.92. Above the sharp lower limit Cmin, 
the eccentricity posterior probability falls off gradually, 
as we discuss in §3. Note that the Cmin we have defined 
here, which assumes we can distinguish between b and 
Peirc (i-e. via some constraint on ingress/egress time), 
is a stronger limit than the minimum eccentricity from 
the constraint that the transit be non-grazing (which we 
discussed in §2.1). 

Returning to the issue of contamination by blending 
(discussed in §2.1), consider a transit with g = 2.5 and 
thus emvn = 0.724. If the transit depth were diluted by 
a factoiQ of 0.9 by an undetected second star in the pho- 
tometric aperture, we would measure g = 0.9^'^2.5 — 
2.435 and infer nearly the same minimum eccentricity of 
Gmin — 0.711. Finally, imagine that some of the con- 
straint on g measured from the light curve came from 
the non-grazing shape of the transit, implying an im- 
pact parameter greater than 1 — Rp/R^,. If the Rp/R^, 

^ This is a worst-case scenario because in fact we could easily 
detect a companion causing such a large dilution. 
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measured from the diluted transit curve were 0.1, the in- 
ferred maximum impact parameter would be 0.9. If the 
true Rp/Ri, is 5% larger, then the maximum impact pa- 
rameter should be 0.895. This translates into a negligible 
effect on the constraint on g. 

In Figure [l] we plot pcirc as a function of uj. Centered 
at w = 90° is a broad range of uj for which pcirc would be 
quite high. For example, for e — 0.9, pdrc would be erro- 
neously high by a factor of 10-100 for -3° < w < 183°, 
over half the possible orientations. Moreover, although 
the periapses of eccentric planets are intrinsically ran- 
domly oriented throughout the galaxy, based on geom- 
etry eccentric planets with uj ~ 90° are more likely to 
transit. For example, from a population of planets with 
e = 0.9(0.95,0.99) and a given orbital separation, we 
would be able to observe 19 (39, 199) times as many 
transiting at periapse as at apoapse. 

Another happy coincidence is that the true stellar den- 
sity is unli kely to be higher tha n the Kepler Input Cat- 
alog (KIC, iBatalhaet al][2010l ) value by a factor of 10. 
The opposite situation is common; a star identified as be- 
ing on the main sequence may actuallybe a low-density 
subgiant or giant (e.g. Mann et al. 2012, Dressing et al. 
2012, in prep). Conversely, there are not many stars with 
the density of lead. Even when precise measurements of 
the stellar density are unavailable, our basic knowledge 
of stellar structure and evolution often allows for con- 
straints on the eccentricity. If there exists a population 
of highly-eccentric Jupiter-sized planets, many of them 
will be identifiable from the light curve alone, i.e. we 
would deduce a large emin- 
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Fig. 1. — The ratio of the circular density to the nominal stellar 
density, Pcirc/p*i required for a circular model to account for the 
transit observables of an eccentric planet. The ratio is plotted as 
a function of the planet's argument of periapse. The solid (dotted, 
dashed) line corresponds to a planet with an eccentricity of 0.95 
(0.9, 0.8). For a large range of periapse angles, one would infer 
a density much larger than the nominal value if one modeled the 
eccentric planet's orbit as circular. 



3. GENERATING AN ECCENTRICITY POSTERIOR 
PROBABILITY DISTRIBUTION 

Through an MCMC exploration — in our case imple- 
mented in the Tra nsit Analysis Package software (TAP, 
Gaza k et al.l [201l[ )-we can not only determine emin but 
impose even tighter constraints on a planet's eccentricity. 
For example, in §2 we stated that a candidate whose cir- 
cular density is consistent with the nominal value could 
have any eccentricity (i.e. for any value of eccentricity, 
there is an w that satisfies g{e,uj) = 1). However, for 
5^1, the eccentricity posterior marginalized over uj will 
be dominated by low eccentricity values, even with a flat 
prior on the eccentricity. For example, if e = 0, any value 
of Lj will satisfy g — 1, whereas only a small range of w 
allow for g = 1 and e > 0.9. Thus, because we expect 
planetary periapses to be distributed isotropically in the 
galaxy, a deduced g = 1 is most likely to truly corre- 
spond to a planet with a low eccentricity. By the same 
argument, the eccentricity posterior corresponding to a 
measured g ^ 1 will peak just above Cmin- 

Of course, the transit proba bility also aff ects the eccen- 
tricity posterior distribution (Burke 2008): an eccentric 
orbit with a periapse pointed towards us (a; = 90°) is ge- 
ometrically more likely to transit than a circular orbit or 
an eccentric orbit whose apoapse is pointed towards us. 
We will discuss how an MCMC exploration automatically 
accounts for the transit probability later in this section. 

3.1. Monte Carlo simulation of expected eccentricity 
and UJ posteriors 

To calibrate our expectations for the output of a more 
sophisticated MCMC parameter exploration, we first 
perform a Monte Carlo simulation to generate predicted 
posterior distributions of e vs. uj via the following steps: 

1. We begin by generating a uniform grid of e and uj, 
equivalent to assuming a uniform prior on each of 
these parameters. 

2. Then we calculate g{e,uj) (Equation |3]) for each 
point (e,aj) on the grid. 



3. We compute 
probng = 



a 1 — e^ 



(9) 



where probjj is the probability of a non-grazing 
transit, for each point {e,uj) (|Winnll20"l0l Equation 
9). We generate a uniform random number be- 
tween and 1 and discard the point if the random 
number is greater than the transit probability. 

4. We calculate the periapse distance -^(1 — e) for 
each grid point and drop the point if the planet's 
periapse would be inside the star (effectively impos- 
ing a physically-motivated maximum eccentricity, 
which is most constraining for small a/R^,). 

5. We downsample to a subset of grid points that fol- 
lows a normal distribution centered on g, with a 
width of cTg/g — 0.1, corresponding to a 30% un- 
certainty in the stellar density. To do this, we cal- 
culate the probability 



prob„ = 



1 



/2^ 



cxp 



[9ie,uj)-gf 



2-.^ 



(10) 
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and discard the point (e,aj) if a uniform random 
number is greater than prob„. 

We pfot the resulting posterior e vs. lo distributions 
in Figure [2] for two a/Ri,, one large and one small, and 
Rp = 0.1. The banana shape of the posterior results 
from the correlation between e and uj (i.e. Equation |4]). 

The posteriors reveal that, rather than being inextri- 
cably entwined with lu, the eccentricities deduced from g 
are well constrained. A pcirc consistent with the nominal 
value (g — 1 with pi, constrained to within 30%) is more 
likely to correspond to a small e (e.g. the probability 
that e < 0.32 is 68.3% for a/R^, = 10 and that e < 0.35 
is 68.3% for a/R^, = 300), while circular densities incon- 
sistent with the normal values {g significantly different 
from unity) have a well-defined minimum e, above which 
the eccentricity posterior falls off gently. For example, for 
g = 2.5 and a/Ri, = 300, the probability that e > 0.69 is 
99%. Furthermore, the eccentricity is likely to be close to 
this minimum eccentricity because the range of possible 
Lu narrows as e — ?> 1. For g — 2.5 and a/Ri, = 300, the 
probability that 0.69 < e < 0.89 is 95%. 

Next we explore how the uncertainty in pi, affects the 
eccentricity posterior, quantifying how "loose" this prior 
constraint can be. In Figure |3l we plot eccentricity 
contours using a/R* = 30 for g = 1 (i.e. consistent 
with circular; bottom) and g = 2.5 (top) for five val- 
ues of ap^ / Pi, assuming a normal distribution and that 
(Tg/g = ^ap^l Pi,. For g = 2.5, the measured eccentric- 
ity is always e — 0.79; it has an uncertainty of tp 06 ^°'' 

ffpj Pi, = 0.01 and Iq 07 ^'^'' '^pJ P* = O-^- Thus the ec- 
centricity remains tightly constrained even for large un- 
certainties in the stellar density. For g = 1, the mea- 
sured eccentricity depends more strongly on the uncer- 

+0.34 " '^ 



tainty: e = 0.031^:^3 fo^ ^p* I P* 



0.01 and e = 0.24to:^8 
for (Jp^/pi, = 0.5. Thus for full, ingress, and egress du- 
rations consistent with circular, a tighter constraint on 
the stellar density allows for a stronger upper limit on 
the eccentricity. However, even for a very poorly con- 
strained Pi,, the posterior reveals that the eccentricity is 
most likely to be small. 

3.2. A Bayesian framework for generating posteriors 

In the Monte Carlo simulation in the previous subsec- 
tion, we used random numbers to select grid points in (e, 
oj) that were consistent with the light curve parameters, 
the prior knowledge of the stellar density, and the transit 
probability. An MCMC fitting routine naturally gener- 
ates such a posterior in eccentricity and uj according to 
the following Bayesian framework. 

Let the model light curve be parametrized by e, uj, Pi,, 
and X, where X represents the additional light curve 
parameters (i.e. orbital period, cos (inclination), radius 
ratio, mid transit-time, limb darkening parameters, and 
noise parameters). Let D represent the light curve data. 
We wish to determine the probability of various e and uj 
conditioned on the data, or pToh{e,uj, pi,, X\D). 

According to Bayes' theorem: 

prob(e,aj, Pi,,X\D) oc prob(I?|e,aj, p*,X)prob(e, w,p*,X) 

(11) 
where the final term represents prior knowledge. 

We assume a uniform prior on all the parameters ex- 
cept Pi,, for which we impose a prior based on the stellar 



parameters and their uncertainties. Therefore, we can 
rewrite the equation as: 

prob(e, UJ, Pi,, X\D) ex prob(Z?|e, uj, Pi,, X)prob(p^) (12) 
Next we marginalize over X and pi, to obtain 

prob(e, aj|Z?) ex / pToh{D\e,uj, Pi,, X)pToh{pi,)dXdpi, 

(13) 
the two-dimensional joint posterior distribution for ec- 
centricity and OJ. The first term under the integral is the 
likelihood of the data given e, uj, pi, and X. Thus a uni- 
form prior on both these quantities naturally accounts 
for the transit probability because 'pToh{D\e,uj, Pi,, X) is 
the transit probability; for certain values of e and uj, the 
observed transit D is more likely to occur. Combina- 
tions of parameters that produce no transits are poor 
models, resulting in a low likelihood of the data. Evalu- 
ation of the likelihood prob(Z?|e, uj, pi,, X) is part of how 
we obtain the parameter posteriors through an MCMC 
exploration, the details of which we describe in the next 
subsection. 
Finally, we can marginalize over uj to obtain 

prob(e|Z?) ex / / proh{D\e,UJ, Pi,,X)proh{pi,)dXdpi,duj 

(14) 
Thus, although stellar density, eccentricity, and uj de- 
pend degenerately on light curve properties (Equation 
[6]), a Bayesian approach to parameter space exploration 
translates a loose prior on the stellar density, prob(p*), 
and uniform priors on the intrinsic planetary values of 
eccentricity and uj, into a tight constraint on the planet's 
eccentricity. 

3.3. Obtaining the eccentricity posterior through an 
MCMC sampling method 

When performing light curve fits with eccentric orbital 
models, it is essential to use an MCMC sampling method, 
or some other algorithm for which the time spent in each 
region of parameter space is proportiona l to the probabil- 
ity. We refer the reader to lBowler et al.l ()2010D (§3) for a 
helpful description of the MCMC method. The MCMC 
method can be used to minimize the x^ (iii the limit of 
uniform priors and Gaussian noise) or to maximize what- 
ever likelihood function is most appropriate given one's 
prior knowledge. In our case, we impose a normal prior 
on Pi, an d account for red noi se using a wavelet-based 
model bv lCarter fc Win"nll2009l Obtaining the eccentric- 
ity posterior through an MCMC sampling method offers 
several advantages: 

1. It naturally allows for marginalization over all val- 
ues of w. For example, in the case of a circular den- 
sity near the nominal value (g ~ 1), the chain will 
naturally spend more time at low eccentricities, for 
which a large range of uj provide a good fit, than at 
high eccentricities, for which only a narrow range 
of UJ provide a good fit. 

2. It reveals and comprehensively explores compli- 
cated parameter posteriors. In particular, some 
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Fig. 2. — Contoured eccentricity vs. lu posteriors from Monte Carlo simulations for representative values of g. The points follow a normal 
distribution centered at the indicated value of g (columns) writh a wridth of 10%, corresponding to a 30% uncertainty in p*. We show the 
posteriors for two values of a/ij* (rows). The black (gray, light gray) contours represent the {68.3, 95, 99}% probability density levels (i.e. 
68% of the posterior is contained within the black contour). Over-plotted as a black-and-white dotted line are histograms illustrating the 
eccentricity posterior probability distribution marginalized over ui. 



of the distributions in Figure [2] and [3] have ba- 
nana shapes, which often cause conventional chi- 
squared minimization algorithms to remain stuck 
in the region of parameter space where they be- 
gan. In contrast, an MCMC exploration will even- 
tually fully sample the p osterior distribution. (See 
iChib fc Greenberglll995l for a pedagogical proof of 
this theorem.) Because of the "banana-shaped" e 
vs. u! posterior for high eccentricities (Figure [2] 
and 131), conventional MCMC algorithms, like TAP, 
require many iterations to converge and fully ex- 
plore parameter space. In our case, we test for 
convergence by plotting e and w each as a func- 
tion of chain link and assess if the exploration 
appears random. We also check to ensure that 
the w posterior is symmetric about w = 90°. 
Asymmetry indicates that the chains have not yet 
converged. We note that the variables ecosw 
and esinoj also have a banana-shaped posterior. 
When feasible, we recommend implementing an 
afRne- invariant code such as emcee that more ef- 
ficiently explores banana-sh aped posteriors (e.g. 
iForeman-Mackev et all I2012D . In gSXH we de- 
scribe how to speed up the fit convergence by using 
g instead of e as a variable while maintaining a uni- 
form prior in e and lo. 

3. It allows us to easily impose priors on certain pa- 
rameters, such as the stellar density. If desired, one 
can impose a prior on the eccentricity. In §4, we 
perform an additional fit for each dataset using a 
Jeffrey'flfl prior on the eccentricity, which is appro- 

^ We use a true Jeffrey's prior prob(e) oc 1/e, which we have 
not normalized because we only consider the ratio of probabilities 
when assessing a jump in an MCMC chain. For the fits in Section 
4, for which e^jin is well above 0, this prior is sufficient. However, 
if e = is a possibility (i.e. for g near 1), the reader may wish 
to use a modified Jefi^rey's prior, prob(e) oc l/(e -|- eo), where eo 
is the noise level. We recommend estimating an upper limit on g 



priate if we wish to avoid assumptions about the 
magnitude of the eccentricity. Here we implement 
the prior through regularization (i.e. as an extra 
term in the jump probability). 

4. It automatically accounts for the transit probabil- 
ity, because jumps to regions of parameter space 
that do not produce a transit are rejected. To ad- 
dress what may be a misconception, we empha- 
size that it is unnecessary — and actually a double 
penalty — to impose transit probability priors on 
the eccentricity or periapse. 

5. It provides uncertainties that are more reliable than 
the estimates based on a simple covariance matrix 
(as obtained from traditional least-squares mini- 
mization) because there is no assumption that the 
uncertainties are normally distributed. The un- 
certainties fully account for complicated parame- 
ter posteriors and correlations. Therefore we can 
be confident in the constraints on pcirc even when 
the ingress and egress are not well-resolved. 

We caution that although this Bayesian framework 
is appropriate for obtaining the posteriors of a single 
planet, selection effects must be carefully considered 
when making inferences about a population. 

3.3.1. Using g as a variable for faster convergence 

Using g (Equation |4]) instead of e as a variable in the 
transit fit model avoids the MCMC having to explore a 
banana-shaped posterior. The g variable allows for faster 
convergence and prevents the chain from getting stuck. 
In order to preserve a uniform prior in e and cj, we must 
impose a prior on g by adding an additional term to the 
likelihood function. Following the Appendix 

from the uncertainty in p^irc a^nd p* and solving Equation (|4]l for 
eo using u = 90°. 
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of iBurke et alj ()2007[ ). the transformation from a uniform prior in e to a prior in g is: 



de 
pToh{g)dg = pToh{e)—dg 
og 



(sir 



, , , 9e sm" uj ( sm" u 
prob(5)=prob(e)— = - 

og 



[^1 + sin a;) ± 2gsina;-\/sin a; — 1 



^/slv? LO-l+g^ig^ 



■ 2 \2 

sm wj 



where we have assumed prob(e) = 1 and for which the -I- corresponds to g > 1 and the — to g < 1. 
Therefore, we add the following term to the log likelihood: 



A/: = In 



sin u) fsin a; — ll 



(l + sin oj) ± 2(7sinaj-\/sin w — 1 + 5^ 



Vsi 



sm Lu 



f [g^ 



■ 2 \2 

sm UJ 



(15) 



(16) 



We demonstrate the use of this variable in §4. We note that in our light curve fits, we use g only to explore 
par ameter space, transfor ming the variable to e in order compute the Keplerian orbit, with no approximations, for 
the lMandelfc Agoll ()2002D light curve model. 

3.4. Obtaining the eccentricity posterior from the circular-fit posterior 

The Monte Carlo exploration in i j3.1l was meant to give us a handle on what the eccentricity and uj posterior 
should look like and how they are affected by uncertainty in pi,. However, one could use a more formal version of 
this exploration to obtain posteriors of eccentric ity and ui directly fro m the posteriors derived from circular fits to 
the light curve, an approach that was adopted bv lKipping et al.l (|2012l ). One could maximize the following likelihood 
for the parameters /o*, e, and u: 



C 



1 [g{e,uj)^pi, - jOcirc]^ _ 1 [p* 

2 a? 2 



-'-A-. measured 



In (prob„ 



(17) 



The first term in the likelihood function demands agreement with the Pcirc derived from the circular fit to 
the light curve. If the pcirc posterior is not normal, one could replace this term with the log of the probability of 
g{e,uj)^Pi, given the pcirc posterior. Note that g{e,ui) can either be computed from the approximation in Equation 
(|3]) or by solving and integrating Kepler's equation to obtain the mean ratio of the transiting planet's velocity to 
its Keplerian velocity over the course of the transit. The second term is the prior on p^^ from the stellar parameters 
independently measured from spectroscopy (or asteroseismology) . The final term is the probability of a non-grazing 
transit (Equation IH]) • If one uses the variable g instead e, one should add Equation ([TB)) to the likelihood. We warn 
that this likelihood function drops constants, so although it can be used to generate parameter posteriors, it should 
not be used 
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to compute the Bayesian evidence quantity. 

In the next section, we demonstrate that this approach 
yields the same eccentricity and w posteriors as directly 
fitting for the eccentricity from the light curve. 

4. DEMONSTRATION: MEASURING THE 
ECCENTRICITIES OF TRANSITING JUPITERS 

To demonstrate that the duration aspect of the pho- 
toeccentric effect allows for precise and accurate mea- 
surements of a transiting planet's eccentricity from the 
light curve alone, we apply the method described in ^ 
to several test cases. In t j4!ll we measure the eccentricity 
of a transiting planet that has a known eccentricity from 
RV measurements. In t ^4.2l we inject a transit into short 
and long cadence Kepler data and compare the resulting 
e and ui posteriors. In §4.3| we measure the eccentricity 
of a Kepler candidate that has only long-cadence data 
available. 

4.1. HD 17156 b: a planet with a large eccentricity 
measured from RVs 

HD 17156 b was discove red by the N e xt 20 Star s 
(N2K) Doppler s urvey ([Fischer et al.l [20051 [20071 ). 
iFischer et al.l ([20071 ) reported that the planet has a large 
orbital eccentricity of e = 0.67 ± 0.08. We identified this 
planet and the releva nt r eferences using e xopl anets . org 
([Wright et all 1201 1[ ). iBarbieri et aLl ([2007D reported 
several partial transits observed by small-telescope 
observers throughou t the Northern H emisphere, and 
IBarbieri et all ([2009D and IWinn et al.l §009) observed 
full transits using high-precision, ground-based photom- 
etry. Here we demonstrate that the planet's eccentricity 
could have been measured from the transit light curve 
data alone. 

We simultan eously fi t thre e light curves (Figure Hj), 
one from Barbicri et al. ( 20 091) and two from I Winn et al.l 
([2009) using TAP (,Gazak etall [20TT[ ). which employs 
an MCMC technique to generate a posterior for each pa- 
rameter of the Mandel and Agol (2002) transit model. 
Time-correlated , "red " noise is accounted for using the 
iCarter fc WinnI ([2009D wavelet-based likelihood function. 
To achieve the 2^ (where N is an integer) data points re- 
quired by the wavelet-based likelihood funct ion without 
excessive zero-padding, we trimmed the first iWinn et al.l 
([2009) light curve from 523 data points to 512 data points 
by removing the last 11 data points in the time series. 
Initially, we fixed the candidate's eccentricity at and fit 
for pcirc with no prior imposed, to see how much it differs 
from the well- measured value of pi,. Then we refitted the 
transit light curves with a normal prior imposed on the 
stellar density, this time allowing the eccentricity to vary. 
In both cases, we treated the limb darkening coefficients 
following the l iterature: we fixed the coefficients for the 
IBarbieri et al.l (2009.1 lig h t curv e and left the coefficients 
free for the [Winn et al.l ([2009() light curves. Following 
IWinn et al.l ([2009() . we also included linear extinction free 
parameters for the two Winn et al. (2009) light curves. 
(The published IBarbieri et al.l 120091 light curve was al- 
ready pre-corrected for extinction.) 

Figure [S| shows posterior distributions from a circular 
fit (top row) and an eccentric fit (botto m row) with a 
prior imposed on the stellar density from lGilliland et al.l 
([20111) . In Figure [6] we compare the posteriors generated 
from a) the eccentric fit to the light curve using 5 as a 



parameter (with a prior imposed to maintain a uniform 
eccentricity prior; Equation ([T6[) to posteriors generated 
using: b) a Jeffrey's prior on the eccentricity, c) e instead 
of (7 as a free parameter (to demonstrate that they are 
equivalent), and d) the likelihood-maximization method 
described in §3.4, using the posterior of pcirc from the 
circular fit. The four sets of posteriors closely resemble 
one another. The computation times were about 1 day 
for the circular fit, about 1 day for the eccentric fit us- 
ing g as a parameter, several days for the eccentric fit 
using e as a parameter, and thirty minutes for the like- 
lihood maximization method of §3.4. Note that the fi- 
nal method requires the best-fitting parameters resulting 
from a circular fit to the light curve, including accurate 
parameter posteriors. We therefore caution against using 
the parameters listed in the Kepler public data releases 
for this purpose because those values are the result of 
a least-squares fit and make the assumption of normally 
distributed parameter uncertainties. However, if one has 
already precomputed circular fits using an MCMC algo- 
rithm that incorporates red noise and limb darkening — as 
we have done for all of the Jupiter-sized KOIs (§5) — the 
final method (§3.4) is advantageous because of the de- 
creased computation time. 
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Fig. 4.— Light curves of HD 17156 from IBarbieri et al.l l|2009h 
(top) and Winn et al. (2009) (middle, bottom). A set of eccentric 
model light curves drawn from the posterior are plotted as solid 
lines. 



Based on the circular fit alone, we would infer 
ff(emin, 7''/2) = 2.0, corresponding to a minimum eccen- 
tricity of Cmin = 0.61. From the eccentric fit, we obtain 
a value of e = 0.7llg'Qg using a uniform prior on the 



eccentricity and e = 0.69 



-1-0.16 
-0.09 



using a Jeffrey's prior. 



Therefore, we could have deduced the eccentricity deter- 
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Fig. 5. — Posterior distributions of e and oj for the HD 17156 
transiting system, with eccentricity fixed at (row 1) and free to 
vary (row 2). Row 1: Left; p* derived from circular fit. The solid 
line marks the nominal value. Right: Posterior distribution for 
eccentricity solving Equation JSJl for uj = (solid line), ui = 45° 
(dashed line), and ui = 90° (dotted line). Row 2: Left: Posterior 
distribution for uj from eccentric fit (i.e. a fit to the light curve 
in which the eccentricity is a free paramete r; solid). Gaussian 
illustrating posterior from lFischer et aLl 120071) RV fit (dotted line). 
Right: Same for eccentricity posterior. 



1.0 



0.5 



0.0 





100 200 
CO (deg) 



0.5 0.6 0.7 0.8 0.9 
e 



Fig. 6. — Left: Posterior distribution for to for a fit to the light 
curve using 3 as a free parameter with a uniform prior on the eccen- 
tricity (sold line) and Jeffrey's prior (dotted line). Posterior distri- 
bution using e instead of oj as a free parameter (dot-dashed line). 
Posterior distribution using method described in §3.4 (dashed line). 
Right: Same as left, for eccentricity posterior. 



mined from 33 RV measurements — e = 0.67 ± 0.08 
(jFischer et al.l I2007D — from these three transit hght 
curves alone. 

The host star has a part icularly well-constrai ned den- 
sity from asteroseismology (jGilliland et al.l[20TTI ). We ar- 
tificially enlarge the error bars on the stellar density from 
1% to 20% and repeat the fitting procedure, obtaining an 



eccentricity of e = 0.70 



+0.14 
-0.09- 



We also repeat the fitting 



procedure with a density derive d from the st e llar p aram- 
eters M^ and i?* determined by I Winn et"al] (|2009f ) from 
isocrone fitting. This "pre-asteroseismology" density has 
an uncertainty of 10% and, r noreover, is abo u t 5% larger 
than the value measured by iGilliland et al.l ()2011[ ). We 
obtain an eccentricity of e = 0.70to;ii- In Figure [7] and 
[51 we plot the resulting posterior distributions, which are 
very similar. Therefore, even with uncertainties and sys- 
tematics in the stellar density, we can measure a transit- 



ing planet's eccentricity to high precision and accuracy. 
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Fig. 7. — Posterior distributions of e and oj for the HD 17156 
transiting system, with three diffe rent priors o n the stellar density: 
the density measured by Gilliland et al.l IJ20111 ) (solid); the density 
measured by Gilliland et al. (2011]) with uncertainties enlarged to 
cfp^lp-k = 0.2, (dashed) an d the density based on the stellar pa- 
rameters from lWinn et aO |(200S ) (dotted). 



4.2. Short vs. long cadence Kepler data 

IKippind ()2010bf ) explored in detail the effects of long 
integration times and binning on transit light curve mea- 
surements, with a particular focus on long-cadence Ke- 
pler data. He demonstrated that by binning a finely- 
sampled model to match the cadence of the data, as TAP 
has implemented, one can fit accurate (though less pre- 
cise than from short cadence data) light curve param- 
eters. Using short and long cadence Kepler data of a 
planet with known parameters (TrES-2-b), he validated 
this approach. 

Here we explore, through a test scenario of an eccentric 
planet injected into short and long Kepler data, whether 
this approach holds (as one would expect) for fitting an 
eccentric orbit and what value short-cadence data adds 
to the constraint on eccentricity. We chose parameters 
for the planet typical of an eccentric Jupiter and main- 
sequence host star: P — 60 days, i — 89.5°, Rp/Ri, = 
0.1, e = 0.8, io — 90°, Af* — Ri, — \, and limb darkening 
parameters /xi = /i2 = 0.3. We considered the situation 
in which long cadence data is available for Q0-Q6 but 
short-cadence is available only for one quarter (or may 
be in the future). We retrieved Q0-Q6 data from the 
Multimission Archive at the Space Telescope Science In- 
stitute (MAST) for Kepler target star KIC 2306756, se- 
lected because it has both long and short cadence data. 
Then we applied the TAP MCMC fitting routine to fit a) 
one short-cadence transit (fixing the period at 60 days) 
that took place in a single segment of short-cadence data 
and b) all seven long-cadence transits. 

As in §4.1, we performed one set of fits fixing the or- 
bit as circular and another set with g and w as free pa- 
rameters, imposing a prior on the stellar density corre- 
sponding to a 20% uncertainty in the stellar density and 
a prior on g from a uniform prior in e and u (Equation 
[T6)) . In both cases, we allowed the limb darkening to be a 
free parameter. We plot the resulting posterior distribu- 
tions of eccentricity and ui in Figure [21 From the circular 
fits, the constraint on pcirc is somewhat stronger from 
the short cadence data (26.3^i[g pq) than from the long 

cadence data (25.9^2'? Po); as IKippind ()2010bD found. 
From the short cadence data, we measure an eccentricity 
of e = 0. 85^0^05 with a uniform prior on the eccentricity 
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1I20T1I) (leffl^ lUiIIiTand et al.l 1I20T1I') with error bars enlarged to 20% (middle); and lWinn etahl l l2009t) (left). 



and e = O.SSto^os with a Jeffrey's prior. From the long 
cadence data, we measure an eccentricity of e = 0.84^005 
with a uniform prior on the eccentricity and e = 0.84_qq4 
with a Jeffrey's prior. Therefore, the long cadence data 
is sufficient to obtain a precise eccentricity measurement. 
In this case, the 20% uncertainty in the stellar den- 
sity dominated over the constraint from the transit light 
curve on pcirc! however, for very well-constrained stellar 
properties, we would expect the greater precision of the 
short cadence data to allow for a tighter constraint on 
the eccentricity (see Figure [3|). 

4.3. KOI 686.01, a moderately eccentric, Jupiter- sized 
Kepler candidate 

KOI 686.01 was ide ntified bv lBorucki et al.l (|2011| ) and 
iBatalha et al.l ()2012[ ) as a 11.1 i?Earth candidate that 
transits its host star every 52.5135651 days. We retrieved 
the Q0-Q6 data from MAST and detre nded the light 
curve using AutoKep f Gazak et al.l 120 ill) . We plot the 
light curves in Figure 1101 

We obtained a spectrum of KOI 686 using the High 
Resolution E chelle Spectrome ter (HIRES) on the Keck 
I Telescope (jVogt et al.lll994D . The spectrum was ob- 
tained with the red cross-disperser and 0'.'86 slit using the 
standard setup of the California Planet Survey (GPS), 
but with the iodine cell out of the light path. The ex- 
tracted spectrum has a median signal-to-noise ratio of 
40 at 5500 1, and a resolution A/AA w 55,000. To 
estimate the stellar temperature, surface gravity, and 
metallicity, we use the SpecMatch code, which searches 
through the CPS's vast library of stellar s pectra for stars 
yfithSpectroscoyy Made Easy {SME: Na\eiit\ fc Piskunovl 
119961 : iValenti fc Fischeill2005[) parameters and finds the 
best matches. The final values are the weighted mean 
of the 10 best matches. We then interpolate these stel- 
lar parameters onto the Padova stellar evolution tracks 
to obtain a stellar mass and radius. We c hecked these 
values using the empirical relationships of iTorres et al.l 
(I2OIOI) . We find p^ = 1.02to:29 Pq (the other stehar 
parameters for this KOI and parameters for other KOI 
will be published as part of another work, Johnson et al. 
2012, in prep). 

We then fit circular and eccentric orbits to the transit 
light curve, as described above, binning the model light 
curves to match the 30-minute cadence of the data. We 
impose a normal prior on the li mb-darkening coefficients 
based on the values from iSind (f201(I) . Figure [11] shows 



posterior distributions from a circular fit (top row) and 
an eccentric fit (bottom row) with a prior imposed on 
the stellar density. We measure the eccentricity to be 

e = 0.62trii 

We cau tion that this cand i date has not yet been 
validated; iMorton fc Johnson! (|2011l ) estimate a false- 
positive probability of 8%. If the candidate is a false 
positive, its orbit (and other properties, such as its ra- 
dius) is likely to be different from that inferred. However, 
we note that if the candidate is a background binary or 
hierarchical triple and is actually larger than a planet, 
the inferred eccentricity would actually be higher (i.e. 
if the candidate is actually larger, it must be moving 
through its ingress and egress even faster), unless KOI 
686 is not the primary and the primary has a higher den- 
sity than KOI 686. Another possibility, if the candidate 
is false positive, is that the assumption of Mp <^ Mi, may 
no longer hold and po (Equation [6|) should be compared 
to Pi, + ^companion rather than pstar to obtain g. How- 
ever, even if /Ocompanion ^ P*, the error in g would be only 
( i)3 = 12.5%. 

iSanterne et al .l (|2012|) recently found a false positive 
rate of 35% for Jupitcr-sized candidates, comprised of 
brown dwarfs, undiluted eclipsing binaries, and diluted 
eclipsing binaries. In the case of diluted eclipsing bi- 
naries, the blend effects that we discus sed in jj2 could 
be larger than we considered. However, iMortonl (I2012D 
notes that most of the false positives that ISanterne et all 
()2012t ) discovered through radial-velocity follow-up al- 
ready exhibited V-shapes or faint secondary eclipses in 
their light curves. In the search for highly eccentric 
Jupitcrs, we recommend a careful inspection of the tran- 
sit light curve for false-positive signatures and, when pos- 
sible, a single spectroscopic observation and adaptive- 
optics imaging to rule out false-positive scenarios. 

If the planetary nature of this object is confirmed, 
it will be one of a number of Jupiter-sized planets 
with orbital periods of 10-100 days and moderate ec- 
centricities, but the first in the Kepler sample with a 
photometrically-measured eccentricity. Many previously 
known, moderately-eccentric planets have orbits inside 
the snow line; their eccentricities are thought to be sig- 
natures of the dynamical process(es) that displaced them 
from their region of formation. 

5. A PLAN FOR DISTILLING HIGHLY-ECCENTRIC 
JUPITERS FROM THE Kepler SAMPLE 
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Fig. 9. — Posterior distributions of e and ui for an injected, artificial transit, with eccentricity fixed at (panel 1) and free to vary (panel 
2-3). The sold curves are from a fit to seven light curves from the long-cadence data and the dotted to a single light curve from the short 
cadence data. Left: p* derived from circular fit. The dashed curve represents the nominal value and its uncertainty. Middle: Posterior 
distribution for tu from eccentric fit (solid line). Right: Eccentricity posterior. 



To test the HEM hypothesis (S12), we are "distill- 
ing" highly-eccentric, Jupiter-sized planets — proto hot 
Jupiters — from the sample of announced Kepler can- 
didates using the publicly released Keple r light curves 
(jBorucki et all \20m iBatalha et all [2012I ). To identify 
planets that must be highly eccentric, we are refitting the 
Kepler light curves of all the Jupiter-sized candidates us- 
ing the TAP. Initially, we fix the candidate's eccentricity 
at 0. We identify candidates whose posteriors for pcirc are 
wildly different than the nominal value p* from the KIC. 
From this subset of objects, we obtain spectra of the host 
stars. We refine the stellar parameters using SpecMatch, 
interpolate them onto the Padova stellar evolution tracks 
to obtain a stellar mass and radius, and check the in- 
ferred M^, a.nd R^ using the empirical relationships of 
[Torres et al.l ()201C1[) . We validate t he candidate using the 
method outlined in iMortoiil ()2012[ ). Finally, we refit the 
transit light curves with a prior imposed on the stellar 
density, this time allowing the eccentricity to vary. This 
process will allow to us easily identify the most unam- 
biguous highly-eccentric hot Jupiters. 

6. DISCUSSION 

Measuring a transiting planet's orbital eccentricity 
was once solely the province of radial-velocity obser- 
vations. Short-period planets were discovered by tran- 
sits and foUowed-up with RVs, which someti mes re- 
veale d a sizable ecce ntricity (e.g. HAT-P- 2b. Bakos et al.l 
120071: CoRoT-lOb, iBonomo et all IMol) . Long-period 
planets — which, based on the RV distribution, are more 
commonly eccentric — were discovered by radial-velocity 
measurements and, on luc ky occasions, fou nd to tran- 
sit (e.g. HD 17156 b, Fi scher et al.l 120071 the planet 
discu ssed in gH as well as HD 806066 b, iNaef et al.l 
1200 ID . But now, from its huge, relatively unbiased tar- 
get sample size of 150,000 stars, Kepler has discovered 
a number of long-period, transiting candidates. Among 
these are likely to be a substantial number of eccentric 
planets (S12), whi c h have enhanced t ransit probabilities 

( 



Kane et al.ll2012D. iMoorhead et al.l (j2011h : iKane et al.l 
20121 ) and lPlavchan et al.l (|2012D have characterized the 



eccentricity distribu tions of these candid ates based on 
Kepler photometry. iKipping et al.l (J2012D are employing 
MAP to measure the eccentricities of planets in systems 
in which multiple planets transits. Here we have demon- 
strated that it is also possible to constrain an individual 
planet's eccentricity from a set of high signal-to-noise 
transits using a Bayesian formalism that employs rela- 
tively loosely-constrained priors on the stellar density. 
The technique we have presented can be applied to any 



transit light curve, as we did in Ji4.1[ for HD 17156 b us- 
in g ground-based photo metry. Comparing this technique 
to lKipping et al.l (j2012D 's MAP, MAP is more model in- 
dependent - requiring no knowledge at all of the stellar 
density ~ but our technique is applicable to single tran- 
siting plan ets, as J upiter-sized Kepler candidates tend 
to be (e.g. iLatham et al. 2011). We are the process of 
fitting the orbits of all Jupiter-sized Kepler candidates, 
which will lead to the following prospects: 

1. For candidates with host star s too faint for RV 
follow-up (65% of candidates in lBorucki et al.|[2011l 
are fainter than Kepler magnitude 14), our tech- 
nique will provide an estimate of the planet's ec- 
centricity. We may also be able to deduce the 
presence of companions from transit timing vari- 
ations, thereby allowing us to search for "smoking 
gun" perturbers that may be responsible for the in- 
ner planet's orbital configuration. In a companion 
paper (Dawson et al. 2012, in prep), we present 
the validation and characterization of a KOI with 
a high, photometrically-measured eccentricity and 
transit timing variations. 

2. For candidates bright enough for follow-up RV 
measurements, the eccentricity and uj posteriors 
from photometric fits allow us to make just a few 
optimally timed radial velocity measurements to 
pinpoint the planet's eccentricity, the mass and 
host-star density, instead of needing to devote pre- 
cious telescope time to sampling the full orbital pe- 
riod. The tight constraints on eccentricity from 
photometry alone can be combined with radial- 
velocity measurements to constrain the candidate's 
orbit — either by fitting both datasets simultane- 
ously or by using the posteriors from the photom- 
etry as priors for fitting a model to the RVs. To 
maximize the information gain, the prior on the 
stellar density should remain in place. This serves 
as an additional motivation for measuring the spec- 
troscopic properties of candidate host stars in the 
Kepler field. 

3. We can also measure the spin-orbit angles of 
the candidates orbiting the brightest stars with 
Rossiter-McLaughlin measurements. Then we can 
compare the distribution of spin-orbit angles of 
those planets we have identified as eccentric with 
the distribution of those we have constrained to be 
most likely circular. 



Photoeccentric Effect 



13 



1.00 



0.98 



0.96 



X 

3 



T3 

CD 
T3 

C 
CD 

^ 0.94 

T3 

T3 

Q> 
_N 

15 



O 



0.92 



0.90 



0.88 



1 — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — r 
KOI-686 




■■» f t fc t t«* » g •»■ 




» ^■ «»■ 



■>• »• — • *■« 



;»«i. 



J I L 



1 



iv. 



J I I L 



;i»\!>:?«s;! 



1 



J I I L 



1 



J I L 



-5 5 

(time - mid transit) [inrs] 

Fig. 10. — Light curves of KOI 686. A set of eccentric model light 
curves drawn from the posterior are plotted as solid lines. The 
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S12 argue that HEM mechanisms for producing 
hot Jupiters should also produce a population of 
highly eccentric (e > 0.9) proto hot Jupiters and 
predict that we should find 3-5 in the Kepler sam- 
ple. Moreover, Kepler^s continuous coverage may 
offer the best prospect for detecting highly eccen- 
tric planets, against which RV surveys are biased 
(| Jones et alJl2006l: 10 'Toole et al.ll2009D . In gH we 
described our process for distilling highly-eccentric 
Jupiters from the Kepler sample. 
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Fig. 11. — Posterior distributions for KOI 686.01 with eccentric- 
ity fixed at (row 1) and free to vary (row 2). Row 1: Left: p* 
derived from circular fit. The solid line marks the nominal value. 
Right: Posterior distribution for eccentricity solving Equation Jsj 
ioT cj = (soUd Une), uj = 45° (dashed line), and uj = 90° (dotted 
line). Row 2: Left: Posterior distribution for uj from eccentric fit 
(solid). Posterior distribution using method from §3.4 (dotted). 
Right: Same as left, for eccentricity posterior. 

The Kepler sample has already revealed a wealth of in- 
formation about th e dynamics and architectures of plan- 
etary systems (e.g. iLissauer et al.l[20Tlt iFabrvckv et al.l 
120121) but primarily for closely-packed systems of low 
mass, multiple-transiting planets. Measuring the eccen- 
tricities of individual, Jupiter-sized planets in the Kepler 
will allow us to investigate a different regime: plane- 
tary systems made up of massive planets that potentially 
underwent violent, mutual gravitational interactions fol- 
lowed by tidal interactions with the host star. 
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